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Ion beam requirements for fast ignition are investigated by numerical simulation taking into 
account new effects such as ion beam divergence not included before. We assume that ions are 
generated by the TNSA scheme in a curved foil placed inside a re-entrant cone and focused on the 
cone apex or beyond. From the focusing point to the compressed core ions propagate with a given 
divergence angle. Ignition energies are obtained for two compressed fuel configurations heated by 
proton and carbon ion beams. The dependence of the ignition energies on the beam divergence 
angle and on the position of the ion beam focusing point have been analyzed. Comparison between 
TNSA and quasi-monoenergetic ions is also shown. 
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I. INTRODUCTION 

Fast ignition (FI) was proposed 20 years ago as an 
alternative to the standard central ignition scheme of in¬ 
ertial fusion targets [l| . It intends to reduce the drive re¬ 
quirements by separating fuel compression and ignition. 
In electron-driven fast ignition (EFI), a fast electron jet 
generated by ultra-high intensity lasers triggers ignition 
of the thermonuclear fuel. Due to the very large diver¬ 
gences and the too high kinetic energies found in EFI ex¬ 
periments and simulations ii, ion-driven fast ignition 
(IFI) has been taking an increasing interest over the last 
years. Ion fast ignition ai offers several advantages 
over EFI, such as generation of collimated beams, well 
known interaction with the plasma and higher flexibility. 
Some exarnples of such a flexibility are the control of ion 
spectra la, [ 3 , the possibility of choosin g th e optimal ion 
species [a, , including deuterium ions [3 , and the use 

of multiple beams for target irradiation [3, mill. IFI 
progress so far can be summarized by the experimental 
achievement of the required ion energies, high conversion 
efficiencies [l^ and ion beam focusing [I^. This last 
achievement is crucial because it points a way forward 
to get beam focusing into the 30-40 /rm spots required 
for IFI. A full review of the current status of IFI can be 
found in Ref. 

Just after the first experimental evidence of proton 
acceleration by the Targ et Normal Sheath Acceleration 
(TNSj^ mechanism [^, its application to FI was pro¬ 
posed Q . This was followed by theoretical studies on the 
TNSA scheme t arget studies [isl - l^ . new irradi¬ 

ation schemes (a Mlii and the use of ions heavier than 
protons [l^, HJ. Other schemes such as the two-pulse 
scheme described in Ref. [IJ have been proposed also. 
This scheme consists of creating a plasma channel fol¬ 
lowed by Deuterium-Tritium (DT) acceleration close to 
the dense core. Hole boring scalings have been studied 
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recently in Ref. [^ . 

Most of the ignition studies carried out so far are based 
on rather ideal target configurations and simplified ion 
interaction models. For instance, these models do not 
take into account cone tip and coronal plasma energy 
deposition nor scattering of the beam ions in the hot 
plasma between the cone tip and the dense core. 

Fully integrated numerical modeling of the IFI scheme, 
from beam generation to fuel ignition, is not possible to¬ 
day with the present computer resources. It requires the 
integration of physical processes with very different spa¬ 
tial and temporal scales. In this paper, we use a sim¬ 
plified model that assumes an ideal initial distribution 
function of ions and computes ion energy deposition by 
solving the Fokker-Planck equation. Specifically, we as¬ 
sume that a highly uniform laser beam, required to im¬ 
prove the ion beam focusing , with irradiances of the 
order of 10^° W/cm^ impinges on a curved foil placed 
inside a cone, generating a proton or carbon ion beam at 
the foil rear surface by means of the TNSA scheme [^ . 
Ions are focused into the inner surface of the cone tip 
or beyond by the electrostatic fields generated near the 
cone walls. From the beam focusing point, ions diverge 
towards the dense fuel. This paper deals with determin¬ 
ing the beam requirements to ignite a target as a function 
of the ion beam divergence taking into account the ion 
interactions with the cone tip and their scattering with 
the background plasma ions. The goal has been to ob¬ 
tain realistic estimates of the ignition energies in the IFI 
scenario. 

This paper is organized as follows. In Section II, the 
computational model used in the simulations is briefly 
outlined. Details of the multidimensional Fokker-Planck 
ion energy deposition model and its differences with the 
standard ion tracking scheme are shown. In Section III, 
ion energy deposition and scattering by the cone tip, en¬ 
ergy deposition of divergent beams in the DT core, igni¬ 
tion energies of proton and carbon ions with maxwellian 
and quasi-monoenergetic energy distributions, and the 
effect of the position of the beam focusing point are pre¬ 
sented. Finally, conclusions and future work are summa- 
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FIG. 1: (a) Sketch of the simulation box used in this work. 
Ions are injected at the inner surface of the cone tip sited 
at z = 5 and propagate towards the dense core. The 
distance from the ion source to the center of the imploded fuel 
is 650 fim. (b) Density profiles along the central line of the 
simulation box. The solid line corresponds to the imploded 
one-shock target of Ref. and the dashed line to an ideal 
super-gaussian density distribution. 


rized in Section IV. 


II. SIMULATION MODEL 

We assume a uniform (constant flux within its cross 
section) ion beam impinging on a compressed DT fuel. 
The beam is generated by the TNSA scheme [l^, which 
yields ions with a maxwellian energy distribution. In 
order to reduce the ion time spread along the stand-off 
distance from the generation foil to the dense fuel, we 
assume that a re-entrant cone is attached to the com¬ 
pressed fuel for generating the ion beam as close as pos¬ 
sible to the DT. Two configurations of the compressed 
fuel at stagnation have been considered. The first one 
was taken from the single-shock one-dimensional (ID) 
target design proposed by Clark&Tabak specifically for 
fast ignition |28|. The simulation box corresponding to 
this case is shown in Fig. [T^. This fuel configuration 
is quite ideal because the distortions induced by the re¬ 
entrant cone in the shell during target implosion and fuel 
compression are not accounted for. It presents, however, 
a large coronal plasma surrounding the core which clearly 
overestimates the plasma densities from the cone apex to 


the DT core when compared with the density profiles ob¬ 
tained in cone-target simulations [2^ . As ignition ener¬ 
gies depend substantially on the areal density of the coro¬ 
nal plasma [2^, the ID imploded target density distribu¬ 
tion shown in Fig. [1^ allow us to obtain an upper limit 
of the ignition energies for the beam conditions studied. 
The imploded fuel areal density is pR = 2.4 gcm“^ and 
we assume an initial temperature of 100 eV in all the 
DT with exception of the small central dip, which has a 
higher temperature in order to keep the pressure balance 
with the surrounding cold and dense fuel. For the ini¬ 
tial density and temperature assumed for the imploded 
core, 450 gcm“^ and 100 eV, respectively, SESAME ta¬ 
bles give a DT pressure of 72 Gbar, which corresponds to 
a DT adiabat relative to the Fermi degenerate pressure 
a = 1.25 {a = p (Mbar)/[2.17 with p in gcm“^). 
A rather low adiabat has been chosen because the target 
was designed for a highly isentropic implosion in order 
to get a quasi-isochoric stagnated fuel distribution [^ . 
Note that the distance between the inner surface of the 
cone tip, where the ion beam is injected, and the core 
center is 145 pm and the distance to the peak density 
around 100 pm , which are similar to those reported by 
Shay et al. [^ . 

The second configuration of the compressed core is the 
super-gaussian density distribution 440exp[—(r/dr)^^], 
where r is the distance to the center, which is placed 
at z = 150 pm, and Sr = 50 pm. This distribution is 
sited on a 10 gcm“^ coronal plasma with a negligible 
areal density, as shown in Fig. [TJd. Hence, a lower limit 
of ignition energies will be obtained in this case. The 
initial temperature is again 100 eV. It is worth noting 
that the super-gaussian density distribution assumed is 
similar to that used in recent electron-driven fast ignition 
studies (sol-l^. 

Diamond-like-carbon (DLC), aluminum, copper and 
gold have been chosen as materials for the cone tip in 
order to check the importance of ion energy losses and 
scattering on the beam energy deposition in the core. In 
realistic target designs, the cone should be protected from 
the x-rays coming from shell implosion by using heavier 
materials, such as lead or gold, at the cone walls. In 
addition, the cone tip should be thick enough to avoid 
rebound shock breakout. 

Calculations have been performed with the radiation- 
hydrodynamics code SARA, that includes Eulerian hy¬ 
drodynamics in cylindrical r-z geometry, flux-limited 
electron conduction, multigroup radiation transport, ion 
energy de pos ition, DT fusion reactions and a-particle 
transport [3J, • Ion energy deposition and a-particle 

transport are computed by the same three-dimensional 
(3D) Eokker-Planck (FP) transport module. Equations 
of state have been taken from SESAME tables . The 
numerical parameters used in the simulations have been 
chosen as follows: cell widths Ar = Az = 1 pm, initial 
time step 3 fs and the number of pseudo-ions injected in 
the simulation box was 2 x 10®. 
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A. Ion energy deposition 

Ion energy deposition is computed by means of the 
kinetic Fokker-Planck (FP) equation. It allows us to in¬ 
clude any ion initial distribution function and the scat¬ 
tering due to ion-ion collisions, neglected in the standard 
ion tracking schemes [2^. The FP equation is solved in 
3D settings by using a Monte Carlo method similar to 
that used in our hybrid code for fast electron transport 
calculations [13 • Here, for simplicity, we consider the 
ID form of the FP collision term, that for suprathermal 
particles reads: 


_ djSilj) rr.d o 

\dt) pp dE dfj, dfi ’ 


( 1 ) 


where tjj = vn is the ion velocity v multiplied by the fast 
ion density n{x,E, ^,t) (’particle angular flux’ in trans¬ 
port theory), S and T are the stopping power and mean 
deflection per unit pathlength, respectively, and /r is the 
cosinus of the polar angle of the ion velocity. The colli¬ 
sion coefficients S and T can be written as 
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where j stands for plasma species, E is the ion kinetic 
energy, rrij, qj, Vj and 0j are the mass, charge, mean ve¬ 
locity and temperature of the background plasma species 
j, m and q are the mass and charge of the fast ion, and 
log A the Coulomb logarithm. 

For ion-ion interactions, we have Xt » I, G{xi) ~ 
1 and E[{xi) ~ 1, while for ion-electron interactions 
in hot plasmas such that v « Ve, one can take the 
limit Xe « I. In this limit, G{xe) ~ 4Xg/(3i/^) and 
E[{xe) « 4xe/(3v^)- With the former assumptions, the 
collision coefficients for an equimolar DT hot plasma can 
be written as: 
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FIG. 2: Range of monoenergetic ions in DT at 400 gem ® as a 
function of the plasma temperature as obtained from Eq.(2). 


ion density. Note that scattering due to ion-ion collisions 
is important for low energies and high electron tempera¬ 
tures 9e , whereas stopping with plasma electrons is dom¬ 
inant for high ion energies and low Of.. 

Ranges of different ion species with typical energies for 
IFI are shown in Fig. [2] as a function of the DT plasma 
temperature 0 {0 = Og = Oi). We have used the stan¬ 
dard stopping power formula for classical plasmas, Eq.(2) 
[sM js^l • This formula predicts range lengthening when 
plasma electron velocities are comparable or higher than 
the fast ion velocity. Range lengthening is important for 
protons with maxwellian energy distributions generated 
away from the DT because it balances the reduction of 
the ion kinetic energy on target due to time spread, keep¬ 
ing the ion range almost constant during the ion pulse 
[201 . It is worth noting that the range lengthening in the 
temperature range O.I — 10 keV is substantially lower 
for ions heavier than protons, as shown in Fig. [21 This 
may be an issue for quasi-mononenergetic ions, as was 
pointed out in previous work [^ [l^l . 

The energy deposition of a perfectly collimated 10 MeV 
proton beam in DT computed by our FP model and by 
the standard ion tracking scheme are compared in Fig. [31 
In both schemes, the energy deposited per unit path- 
length decreases with the ion penetration because the 
electron stopping power is dominant and scales, roughly, 
as for 0i = 0e = 10 keV. The FP energy deposition 
profiles show beam scattering, straggling and blooming 
near the end of the range, the absence of the Bragg’s peak 
and a slight range shortening due to the proton scattering 
with the background ions. In the case of ion tracking, the 
energy deposition peaks at the end of the range, showing 
the well known Bragg’s peak, which may lead to under¬ 
estimated ignition energies. Note that the Bragg’s peak 
disappears in the FP energy deposition profile due to 
proton scattering. For temperatures lower than 10 keV, 
scattering is less important and the differences between 
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FIG. 3: Comparison of the energy deposition of a perfectly 
collimated 10 MeV proton beam in a uniform DT plasma at 
140 gcm“^ and 10 keV (a) standard linear energy deposition 
model (b) Fokker-Planck (FP) model. Protons are injected 
at 2 = 5 pm and pass through a 25 pm thickness DLC layer 
at 3 gcm“® mimicking the cone tip. The energy deposition 
per unit path length obtaned by both models are compared 
in the inset. 


the FP and tracking energy deposition profiles are lower, 
but still significant. 

To assess the importance of the FP energy deposition 
in FI targets, we have compared the ignition energies of 
the DT isochoric sphere initially at 400 gcm“^ and 100 eV 
heated by a uniform cylindrical beam. Atzeni et al. [T^ 
reported for this target a minimum ignition energy of 8.5 
kJ. Our calculations show the same ignition energy when 
the ion scattering is off, while it increases slightly to 8.75 
kJ for the FP energy deposition model. Thus, despite ion 
scattering effects seem not to be important for ion energy 
deposition, they are not negligible and should be taken 
into account in IFI calculations. Heavier ions, such as 
carbon, are less sensitive to scattering, being the ignition 
energies almost the same with and without scattering. 


B. Ion pulse on target 


Beam power on target for proton and carbon ions with 
maxwellian energy distributions are plotted in Fig. 3] 
The beam temperatures are Tp = 7 MeV and Tc = 200 
MeV, respectively. The pulse of quasi-monoenergetic car¬ 
bon ions with a mean energy {]^ = 650 MeV and an en¬ 
ergy spread of 6E/E = 0.125 Q is also shown for com¬ 
parison. All beams have an energy of 10 kJ and are 
generated at a distance to the simulation box d = 500 
/rm. The power on target p{t) of an ion beam generated 
instantaneously in a foil placed at a distance d is given 
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FIG. 4: Beam power at the left surface of the simulation box 
as a function of time for different ion beams with an energy 
of 10 kJ. The distance from the ion-generation foil to the left 
surface of the simulation box is d = 500 pm in all cases, r 
stands for the ion pulse duration on the foil. 


by the formula [^ : 

where to = d{m/2Ti,)^^^, m and Tf, are the beam ion mass 
and temperature, respectively, Ef, the beam energy and 
d the distance between the ion-generation foil and the 
simulation box. From Eq.(9), the beam power on target 
can be obtained for a pulse of duration r. Assuming that 
ions are generated at the source foil at a constant rate 
over the pulse duration r, the power on target P(t) is 
given by 


P{t) = [ 9{t')p{t - (10) 

with t > T and g{t) = Ijr. This formula has been used 
to calculate the pulses shown in Fig. H) We assume an ion 
pulse duration of 1 ps throughout this paper, with excep¬ 
tion of the quasi-monoenergetic beams, for which longer 
pulses, e.g. 3 ps in Fig. [d] are assumed to avoid very 
high peak powers. The feasibility of ps pulses for quasi- 
monoenergetic ion schemes such as BOA (break-out after 
burner) or RPA (radiation pressure acceleration) remains 
to be demonstrated Q. 


C. Beam focusing 

In addition to the ’conventional’ techniques of ion 
beam focusing such as ballistic transport [H, , mag¬ 
netic lenses and self-generated fields in hol¬ 

low microcylinders illuminated by intense sub-picosecond 
laser pulses [d^, laser-driven proton beam focusing has 
been demonstrated experimentally over the last years. 
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Kar et al. showed beam focusing by using rectan¬ 
gular or cylindrical hollow lens attached to a foil target. 
Offermann et al. found theoretically and experimen¬ 
tally that ion divergence depends on the thermal expan¬ 
sion of the co-moving hot electrons, resulting in a hyper¬ 
bolic ion beam envelope. Using these results, Bartal et 
al. demonstrated experimentally an enhanced focusing of 
TNSA protons in cone targets, predicting spot diameters 
about 20 fim for IFI conditions , well under the 40 ^m 
spots required Q. The focusing mechanism reported by 
Bartal et al. is based on the onset of an electron sheath 
near the cone walls. This sheath generates focusing elec¬ 
trostatic fields that avoids the expansion of co-moving 
electrons. As these fields are generated near the cone 
walls, they affect mainly to the low energy ions, keeping 
the initial divergence of the higher energy ions generated 
closer to the beam axis. As a result, even in the case 
of a very good focusing by the cone walls, the beam will 
follow a hyperbolic path when traveling towards the core 
[13 . m, [43 . Qiao and coworkers have pointed out by im¬ 
plicit PIC simulations that the setting up of the electron 
sheath at the cone walls may reduce substantially the 
laser-to-proton conversion efficiency for the long pulses 
required in IFI . This can be mitigated by special de¬ 
signs of cone walls including insulator materials in order 
to reduce the electron flow between target and cone. 

Recent experiments carried out at laser intensities of 
Ri 5 X 10^° Wcm“^ relevant for IFI have evidenced the 
generation of highly collimated proton beams, i.e. with 
divergence full-angles lower than 10°, in planar foils il¬ 
luminated by high-contrast (r; 10“^^) laser beams 
This result together with those reported in Refs. [T3.lT7lj. 
shows a clear path to achieve the focusing requirements 
of IFI if uniform and high contrast laser beams impinging 
on a curved target placed inside a cone are used to ignite 
the fuel. 

As a precise determination of the beam path requires 
detailed PIC calculations, we have assumed that ions are 
focused to the cone apex or beyond and diverge from 
there to the compressed core with a given angle, which 
is taken as a parameter. The reference case calculations 
shown in the next Section assume that the ion beam is 
focused into a 30 fj,ia spot diameter at the inner surface of 
the cone tip. From this spot, the ion beam diverges with 
an opening half-angle /3 as a parameter ranging from 0 
to 20°. For a given (3, each ion is injected in the simula¬ 
tion box with a randomly selected initial angle between 0 
and p. This is an approximation of the curved ion beam 
trajectories found theoretically and experimentally [4^, 
but allows us to assess the divergence effects on beam 
energy requirements. As the assumption of beam focus¬ 
ing at the inner surface of the cone tip may appear to 
be quite restrictive, we have considered also the cases of 
ions focused on the outer surface of the cone tip and on 
the coronal plasma near the dense DT core. In these 
two cases, the divergence effects are much lower than in 
the reference case due to the shorter stand-off distance 
between the focusing point and the core. 



Ion temperature [MeV/u] 

FIG. 5: Ignition energies versus ion temperature for proton 
and carbon ion beams impinging perpendicularly on the tar¬ 
get shown in Fig. [T] The beams diameter is 30 fim. Solid lines 
correspond to the ID imploded target density distribution and 
dashed lines to the super-gaussian density distribution. 


III. RESULTS 
A. Optimal ion temperature 

We start by determining the ion temperatures that 
maximize the beam-target coupling efficiency or, equiv¬ 
alently, minimize the beam energies required to ignite 
the target shown in Fig. [TJ The results are depicted in 
Fig. [ 5 ] for proton and carbon ions with maxwellian en¬ 
ergy distributions. We assume that TNSA carbon ions 
are generated by using thin, micrometer scale DLC foils 
with the rear surface free of protons. The ignition ener¬ 
gies have been obtained as the minimum ion beam energy 
for which the thermonuclear fusion power has an expo¬ 
nential or higher growth sustained in time. 

The shape of the curves shown in Fig. [5] can be ex¬ 
plained as follows. For low temperatures, ions do not 
penetrate enough in the target to deposit a significant 
fraction of their energy into the dense core, raising the 
ignition energies. For high temperatures, ions penetrate 
into the core more than the optimal areal density of 1.2 
gcm“^ 0 raising again the ignition energies due to the 
increase of the DT mass heated by the beam. 

Ignition energies depend on the target density distribu¬ 
tion and, in particular, on the coronal plasma surround¬ 
ing the core. For the ID imploded density profile with 
the long coronal plasma shown in Fig. [Ho, the ignition en¬ 
ergies are minimized for the ’optimal’ ion temperatures 
of Tp = 8 MeV and Tc = 200 MeV for proton and carbon 
ion beams, respectively. For the super-gaussian density 
profile almost without coronal plasma shown in Fig. [TId, 
the optimal ion temperatures are Tp = 6 MeV and Tc 
— 150 MeV, respectively, substantially lower than those 
found for the imploded target density profile. Note also 
that even for the optimal ion temperatures, the ignition 
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beam diameter [^m] 

FIG. 6: Ignition energies of the target shown in Fig. [T]heated 
by a perfectly collimated maxwellian proton beam with tem¬ 
perature Tp = 7 MeV. Protons are injected at the inner sur¬ 
face of the DLC cone tip. Red curve labeled by squares corre¬ 
sponds to the ID imploded target density distribution and the 
blue curve labeled by diamonds to the super-gaussian density 
distribution shown in Fig. [T)d. 

energies are, roughly, 30% higher for the imploded target 
due to the energy deposition in the coronal plasma. Even 
in this case, assuming a laser-to-ion conversion efficiency 
of 15% similar to that obtained recently in proton accel¬ 
eration experiments [l^ and close to the 12% found by 
Snavely et al. [l5j| . the laser beam energy requirements 
are around 100 kJ. 

Independently of the DT density distribution, it is re¬ 
markable that the minimum ignition energies are similar 
for proton and heavier ions such as carbon, in agreement 
with our results for quasi-monoenergetic beams reported 
recently Q. On the contrary, the optimal ignition ener¬ 
gies increase with the ion atomic number Z scaling, ap¬ 
proximately, as (obtained for proton, carbon and 

other ion beams not shown), which is weaker than the 
dependence found for quasi-monoenergetc ions Q. 

In order to simplify the analysis of the next sections, 
we assume Tp = 7 MeV and Tc = 200 MeV as ’opti¬ 
mal’ temperatures for the two density profiles depicted 
in Fig. [Hd. Note that these ’optimal’ temperatures are 
in the plateau region where the ignition energies depend 
only weakly on the ion temperature for both density dis¬ 
tributions, as shown in Fig. [S] 


B. Optimal beam radius 

In order to determine the ion beam focusing require¬ 
ments, we have studied the dependence of the ignition 
energies on the beam diameter. A perfectly focused pro¬ 
ton beam with the ’optimal’ temperature Tp = 7 MeV 
impinging on the target sketched in Fig. [T] has been cho¬ 
sen as model problem. The results are shown in Fig. [51 


For beams with small diameter (< 10 - 15 /rm, depending 
on the density distribution), the energy density deposited 
is very high, heating the DT fuel up to high temperatures 
with the subsequent range lengthening of the beam ions. 
This, together with the strong hydrodynamic response of 
the DT plasma to the high energy deposition may result 
in ions passing through the dense core and even escap¬ 
ing by the target rear surface. As a result, the ignition 
energies rise up to high values. 

For beams with diameter greater than 35 - 40 ^m, de¬ 
pending on the density distribution, the ignition energies 
rise again due to the increase of the DT volume heated by 
the beam. In this case, the ignition energies are sensitive 
to the ion energy losses in the coronal plasma, result¬ 
ing in a substantial difference between the imploded and 
super-gaussian fuel density profiles, as shown in Fig. [51 
Thus, for the simulation conditions analyzed here, beam 
diameters on DT core should be in the range between 
15 and 35 /rm in order to get ignition energies close to 
the minimum value. This range defines the focusing re¬ 
quirements for IFI and is quite similar to that obtained 
in previous works for quasi-monoenergetic ions Q. We 
have assumed a beam diameter of 30 /rm in the analysis 
carried out in the next sections. 


C. Energy deposition 

As it was mentioned in the introduction, ions have the 
advantage of their classical interaction via Coulomb scat¬ 
tering with the background plasma. Unlike EFI, collec¬ 
tive interactions due to self-generated fields do not play a 
significant role in IFI. This can be explained by the lower 
current densities found in IFI, which are about two orders 
of magnitude lower than in EFI for the same particle den¬ 
sity. In addition, for ignition-scale targets, temperatures 
at the cone tip are high enough to mitigate strongly the 
growth of resistive helds. Hence, collective effects such as 
anomalous ion energy deposition, beam self-focusing and 
filamentation of co-moving electrons may not be impor¬ 
tant in IFI. Nevertheless, detailed investigations about 
those effects in IFI conditions are still open. A first step 
in this direction has been given by Qiao et al. [i^ . 

In addition to the initial ion beam divergence at 
the ion-generation foil, ions are scattered when passing 
through the cone tip, increasing its divergence |50|. The 
effect of the cone tip material on the proton beam energy 
deposition is analyzed in Fig. [71 where the energy den¬ 
sity deposited by the target of Fig. [T^ heated by a pro¬ 
ton beam with the ’optimal’ temperature Tp = 7 MeV 
is shown. The increase of the beam diameter near the 
dense core seen in Fig. [7^ evidences the proton scatter¬ 
ing by the cone tip. As the mean deflection coefficient 
T is proportional to the ion charge squared, Eq.(3), fast 
ion scattering is much lower for light materials, such as 
DLC or aluminum, as shown in Fig. |71 d. On the con¬ 
trary, for gold cones, scattering is important resulting 
in a substantial increase of the beam diameter near the 
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FIG. 7: (a) Energy density isocontours of a proton beam with 
a temperature Tp = 7 MeV impinging perpendicularly on the 
target shown in Fig.[T^. The solid gold cone tip thickness is 25 
/im. The white circle shows the initial position of the density 
isocontour p = 200 gcm“®. (b) Lineouts of the energy den¬ 
sity along the dashed line placed at z = 108 pm for different 
cone tip materials. The lineout corresponding to the standard 
proton linear tracking scheme is also shown for comparison. 


dense core. For instance, in the case with a gold cone tip 
of Fig. [7^, the beam radius increases by around 8 pui, 
which is of the same order than that obtained by Key 
et al. [13 by means of the Moliere’s multiple scatter¬ 
ing theory (which provides the scattering angle in solid 
matter, being a lower limit in plasmas [50j). Moliere’s 
theory gives a 3.7° scattering angle for monoenergetic 
protons with energy (3/2)Tp = 10.5 MeV when travers¬ 
ing a 25 pm solid gold foil. The corresponding increase 
of the beam radius at the DT core sited at a distance of 
75 pm from the outer surface of the cone tip is 4.8 pm, 
which is almost a half of the 8 pm shown in Fig. [7 )d. This 
can be explained by the maxwellian, not mononenergetic, 
energy distribution of protons and also by the more ef¬ 
fective scattering by plasma ions than by solid matter. 
Enhancement of beam radii at the DT core results in 
higher ignition energies Eig (see Fig. HI). For instance, 
the ignition energy of the target sketched in Fig. [T] is Eig 
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FIG. 8: Maps of the energy density deposited in the target 
shown in Fig. [T^ by (a,b) protons with Tp = 7 MeV; (c,d) 
carbon ions with Tc = 200 MeV and (e,f) 650 MeV carbon 
ions with an energy spread SE/E = 12.5%. Plots (a), (c) 
and (e) correspond to perfectly collimated ion beams, whereas 
plots (b), (d) and (f) correspond to divergent beams with an 
opening half-angle P = 20°. The white curves show the initial 
position of the density isocontour p = 200 gcm“®. Ignition 
energies Eig for each case are shown. 


= 16.5kJ for DEC and aluminum cone tips, whereas it 
raises to 17.5 kJ and 21.5 kJ for copper and gold cones, 
respectively. Thus, we can conclude that the use of light 
materials at the cone tip is more suitable for IFI, sim¬ 
ilarly to EFI [^. A DEC cone tip is assumed in the 
reminder of this paper. 

The energy deposition by proton and carbon ions with 
the optimal temperatures are compared in Fig. |S] for per¬ 
fectly collimated and divergent beams. Energy deposi¬ 
tion by quasi-monoenergetic carbon ions has been in¬ 
cluded also for comparison. We analyze first the per¬ 
fectly collimated beams, Figs. c and e. Protons with 
a maxwellian energy distribution show a more localized 
energy deposition than carbon ions with the same dis¬ 
tribution. This can be explained by the better balance 
between range lengthening (which is more pronounced for 
protons than for carbon ions) and energy spread found for 
proton beams. This effect is even more evident by com¬ 
paring the energy deposition by maxwellian and quasi- 
monoenergetic carbon ions, which have a very small en¬ 
ergy spread and, thus, their energy is deposited in a vol- 

























ume determined mainly by the range lengthening effect. 
The higher energy deposition by the maxwellian ions in 
the coronal plasma is also remarkable. 

Note in Fig. [8] that the beams have a penetration into 
the dense core around 1.5 gcm“^, and slightly higher for 
the quasi-monoenergetic ions. These values are higher 
than the optimal value of 1.2 gem’ BlEIl due to the 
energy deposition in the coronal plasma. The coupling 
efficiencies, defined as the energy deposited in the DT at 
densities higher than 200 gcm“^, for the perfectly col¬ 
limated beams are 0.58 for the maxwellian proton and 
carbon beams, and 0.70 for the quasi-monoenergetic car¬ 
bon ions, similar to those reported in Refs. [i,|2|]. Cou¬ 
pling efficiencies rise substantially for the super-gaussian 
DT density distribution of Fig. [T]3 due to the lack of 
the coronal plasma, reaching the values of 0.86 for the 
maxwellian proton and carbon ions, and 0.87 for the 
quasi-monoenergetic carbon ions. 

The energy deposition maps of beams with a diver¬ 
gence half-angle jS = 20° are shown in Figs. |B}d, d and 
f. Note the large hot spot elongation due to the high 
beam divergence assumed, which distorts substantially 
the standard rectangular hot spot shape. As a result, the 
energy deposition volumes are higher than those found 
for collimated beams with the subsequent lower coupling 
efficiencies and higher ignition energies. For instance, 
the coupling efficiencies are 0.42 for maxwellian proton 
and carbon ions, increasing up to 0.59 for the quasi- 
monoenergetic carbon ions. Again, these ions penetrate 
deeper in the core, heating a much larger volume than 
the maxwellian ions. 


D. Minimum ignition energies 

The ignition energies for different beams as a func¬ 
tion of the ion divergence half-angle are shown in Fig. [HI 
These energies have been obtained as described in Sec¬ 
tion UFA. As expected, the ignition energies show a rapid 
increase with the beam divergence angle. The small dif¬ 
ferences between maxwellian proton and carbon ions can 
be explained by the lower range lengthening of carbon 
ions, which does not balance completely their energy 
spread along the source-target distance. It is relevant 
here to recall that ignition energies are almost indepen¬ 
dent of the beam ion atomic number, provided that ions 
have the ’optimal’ mean kinetic energy, as can be seen 
in Fig. [5] The much lower ignition energies found for 
quasi-monoenergetic ions show the strong dependence of 
ignition energies on ion spectrum and the better coupling 
of quasi-monoenergetic ions with the plasma. 

For maxwellian ions, the solid line curve labeled with 
solid squares corresponds to the ignition energies for the 
ID imploded target density profile, while the dotted line 
curve labeled by empty squares corresponds to the super- 
gaussian density profile of Fig. [TJd. Both curves show the 
limit cases of high and low coronal plasma areal density, 
respectively. Actually, the ignition energies of a realistic 
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FIG. 9: Ignition energies of the target shown in Fig. [T] 
heated by maxwellian proton and carbon ions with tem¬ 
peratures Tp = 7 MeV and Tc = 200 MeV, respectively. 
The ignition energies with quasi-monoenergetic carbon ions 
are also shown for comparison. Solid lines correspond to 
the ID imploded density profile shown in Fig. [Tb, whereas 
dashed lines correspond to the super-gaussian density profile. 
The curves labeled by squares show the ignition energies for 
maxwellian protons, the line labeled by diamonds corresponds 
to maxwellian carbon ions, and the lines labeled by triangles 
correspond to quasi-monoenergetic carbon ions. 

case would lie in between these curves, within the grey 
shaded area of Fig. [5] If we assume a laser-to-ion con¬ 
version efficiency of 15% [l^ and an ignitor laser beam 
energy of 110 kJ, Fig. [9] shows that target ignition would 
be possible with maxwellian ions if the beam is perfectly 
collimated or has a divergence half-angle lower than 8° 
for high and low coronal plasma areal densities, respec¬ 
tively. 

For quasi-monoenergetic ions, the ignition energies are 
substantially lower than those obtained for maxwellian 
ions. They are within the bottom shaded area of Fig. [HI 
Note that for the case of low areal density, the kinetic 
energy of carbon ions has been reduced from 650 to 400 
MeV in order to have the optimal beam-target coupling. 

For the divergence half-angles of 5° found in planar 
foils experiments (4^ and assuming a laser-to-ion con¬ 
version efficiency of 15%, our results show that it would 
be sufheient with an igniting laser beam energy between 
130 and 90 kJ for maxwellian ions, and from 75 to 60 kJ 
for quasi-monoenergetic ions, depending on the coronal 
plasma areal density. 

E. Beam focusing beyond the cone tip 

In the cases analyzed in Figs. 5 to 9, we have assumed 
that the beam is focused by the cone walls on the inner 
surface of the cone tip. In this case, the stand-off distance 
between the focusing point and the DT core is maximum 
and the beam has more room to spread out. As there is 
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FIG. 10: Ignition energies of the target shown in Fig. [T^ 
heated by maxwellian protons with temperature Tp = 7 MeV. 
Curve labeled by squares shows the ignition energies for pro¬ 
tons injected at the inner surface of the cone tip. The curve 
labeled by diamonds corresponds to protons injected at the 
outer surface of the cone tip and the curve labeled by trian¬ 
gles to protons injected at z = 65 /rm. In all cases the beam 
spot diameter is 30 /rm. 



experimental and simulation evidence of beam focusing 
beyond the cone tip [13, , it is worthwhile extending 

the study carried out in the last Sections to situations 
in which the beam is focused closer to the dense core. 
In this case, the DT heating by divergent beams is more 
efficient due to the lower beam radius at the core. 

In order to analyze the effect of the beam focusing 
point on the ignition energies, two additional cases have 
been considered: (i) a proton beam focused on the outer 
surface of the cone tip and (ii) a proton beam focused 
40 /im away from the dense core (z = 65 /im). In both 
cases, protons are focused on a 30 p,m diameter spot and 
have a maxwellian energy distribution with Tp = 7 MeV. 
The ignition energies obtained for the ID imploded target 
density profile are plotted in Fig. [TUI Note the large 
reduction of the ignition energies when the proton beam 
is focused closer to the DT core. For instance, if the 
beam is focused 40 /im away from the core, the ignition 
energies Eig are much lower than in the reference case 
and depend weakly on the divergence half-angle. Even if 
the beam is focused on the outer rather than on the inner 
surface of the cone tip there is a substantial reduction of 
Eig, specially for large divergence half-angles. 


IV. CONCLUSIONS 

Our simulations have evidenced that ion beam ignition 
energies have been underestimated so far by using sim¬ 
plified energy deposition models and target designs. Fast 
ion spectrum, focusing spot diameter and position, beam 
divergence and ion scattering with the cone tip are key 


parameters that determine ignition energies. Taking into 
account all these parameters requires a kinetic treatment 
of ion energy deposition such as the linear Fokker-Planck 
model used in this work. Thus, our results generalize 
those obtained with the standard linear particle tracking 
method used so far. 

Ignition energies of compressed targets have been eval¬ 
uated in the limits of high and low coronal plasma areal 
densities. The goal has been to obtain an upper and 
lower limits of beam energy requirements for IFF The 
results of our study can be summarized as follows: 

i) The ion beam diameter on the dense core should not 
exceed around 40 /rm, which is even hig her than the fo¬ 
cal spots inferred from experiments |l4| and obtained in 
large scale PIC simulations of ion acceleration in the IFI 
scenario [131 . Beam focusing is one of the most challeng¬ 
ing issues of IFI, with high chances of success in future 
IFI experiments. 

ii) The cone tip material should be as light as possible 
in order to reduce ion scattering. Intermediate materials 
such as copper are acceptable because they give rise to a 
modest increase of the ignition energies when compared 
with light materials. 

hi) As expected, ignition energies depend strongly on 
the beam divergence. Our calculations show that for the 
beam divergence half-angle of 5° measured in the recent 
experiments, the ion beam ignition energies of the tar¬ 
gets considered here are from 13 to 20 kJ, depending on 
the density distribution of the imploded fuel. These en¬ 
ergies are substantially higher than the standard ignition 
energies of 10 kJ considered so far. 

iv) Assuming the laser-to-proton conversion efficiency 
of 15% found in recent experiments [l^, the laser beam 
energies required for ignition are in between 90 and 130 
kJ for small and large coronal plasma areal densities, re¬ 
spectively. In principle, these energies could be affordable 
with the present laser technology. 

v) The position of the ion beam focusing point is cru¬ 
cial to reduce the ignition energies. If it is beyond the 
cone tip, the ignition energies would be reduced substan¬ 
tially and would be less sensitive to the ion divergence. 

vi) If beam focusing were demonstrated for advanced 
schemes, such as the break-out afterburner scheme [1], 
quasi-monoenergetic ions instead of TNSA ions should 
be used due to their better coupling with the compressed 
fuel. 

The challenge over the forthcoming years will be to 
perform experiments relevant for a full characterization 
of ion generation, focusing and interaction with imploded 
targets in conditions close to those found in IFI. The 
success of IFI will depend on the possibility of generat¬ 
ing ions by means of high contrast, highly uniform laser 
beams and focusing the ion beam into 20 - 30 /rm spots 
as close as possible to the imploded DT core. 

Future studies will include PIC simulations of ion ac¬ 
celeration and focusing to get a full characterization of 
the ion source. This will be followed by integrated sim¬ 
ulations similar to those presented here, but with a self- 
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